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Abstract 

We investigate the performance of Grover's quantum search algorithm on a register which is 
subject to loss of the particles that carry the qubit information. Under the assumption that the 
basic steps of the algorithm are applied correctly on the correspondingly shrinking register, we 
show that the algorithm converges to mixed states with 50% overlap with the target state in 
the bit positions still present. As an alternative to error correction, we present a procedure that 
combines the outcome of different trials of the algorithm to determine the solution to the full 
search problem. The procedure may be relevant for experiments where the algorithm is adapted as 
the loss of particles is registered, and for experiments with Rydberg blockade interactions among 
neutral atoms, where monitoring of the atom losses is not even necessary. 

PACS numbers: 03.67Lx, 03.67.Pp, 03.65.Yz 



I. INTRODUCTION 



The quantum search protocol over an unstructured database of N = 2 n elements (n is the 
number of qubits) identifies the correct n-bit string that matches a given target state. The 
search algorithm originally proposed by Grover pQ is an iterative algorithm, where each step 
in the algorithm involves two operations: A change of sign of the target state amplitude and 
an inversion of all quantum state amplitudes about their mean in the computational basis 
[TJ |2]. The ideal evolution during a single iteration of the algorithm can be described by 
a rotation in the two-dimensional subspace spanned by the target state and the symmetric 
superposition of all logical states. This explains the functioning of the algorithm and allows 
exact evaluation of the unitary dynamics. After j\/N iteration steps, the algorithm has 
rotated an initially symmetric superposition of all computational basis states of the qubit 
register into a state very close to the target state. The advantage of Grover's algorithm 
(GA) over a classical search algorithm is obtained only when the size of the search space 
is large (N ^> 1). The larger the database, the longer, however, is the time (~ y/N) for 
finding the target state, and the more important are the effects of error processes relevant 
to the physical system [3H7]. While error processes that preserve the two-level picture of 
the quantum search have been studied by many authors [8j [9] , single qubit dephasing and 
dissipation leak population out of this subspace and one has to resort to numerical methods 
for estimating the success probability of the algorithm [TUHT2"] . 

An important error process in some physical schemes for quantum computing is the 
physical loss of the particles that carry the qubit information. Such errors are relevant for 
neutral atom-based quantum computing protocols, where losses due to collisions with the 
background gas and due to the heating by laser excitation are common at long time scales 
[13] [H]. Schemes to identify and correct errors due to losses have received some atten- 
tion. The fact that we may know precisely which qubit is lost and has to be replaced leads 
to simpler schemes and quite high error thresholds for fault tolerant quantum computing 
compared to the case of general qubit errors [T5HT8] . In this work we shall provide a full 
theoretical analysis of Grover's algorithm in the presence of random qubit losses during the 
search operation. In our analysis we shall make no attempts to detect or correct particle 
losses. Instead, we will assume an implementation of the algorithm where the sign changes 
are implemented perfectly on the relevant Hilbert space of the remaining qubits, i.e., a sign 
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FIG. 1: (a). Illustration of physical loss of qubits from a quantum register, (b). Schematic repre- 
sentation of the block-diagonal density matrix, with population leaking from one two-dimensional 
subspace to the other caused by the loss of qubits. 

change is made on the state component with correct bit values in all the positions of the 
remaining qubits and the amplitude inversion about the mean refers to the state amplitudes 
in the computational tensor product basis of the remaining qubits. As we will show, the 
Grover steps proceed while the system passes through a sequence of two-dimensional sub- 
spaces of the shrinking tensor product spaces, and eventually we end up with a candidate 
solution for a subset of the target bit values. We suggest an effective method by which 
one can reconstruct the full target state after a modest number of independent trials have 
provided candidate solutions for all bit values. 

The paper is structured as follows. In Sec. II, we shall describe the effect of qubit loss 
as a transition between two-dimensional subspaces of the tensor product spaces of different 
numbers of qubits. In Sec. Ill, we describe the Grover steps by an effective Hamiltonian. We 
then derive a continuous time master equation which describes the combined effect of the 
algorithm and the random loss process and we analyze the resulting dynamics. In Sec. IV, 
we describe a method to efficiently obtain the solution of the search problem from different 
subsets of surviving qubits in independent trials of the algorithm. In Sec. V, we describe a 
particular physical implementation of quantum computing using neutral atoms and Rydberg 
blockade interactions, for which the present analysis applies. Sec. VI concludes the article. 
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II. GROVER ALGORITHM IN THE PRESENCE OF LOSSES 



In this section, we show explicitly that the Grover algorithm, both in the absence and 
in the presence of losses can be described by quantum dynamics within two-dimensional 
subspaces. Each loss event merely shifts the dynamics from one 2D space to another. 
Henceforth, the total evolution will be described by a concatenated evolution among these 
subspaces. Since the remaining qubits are entangled with the lost ones, they are properly 
described by a mixed state with a corresponding reduced density matrix. We thus loose 
both the information carried by the lost qubits and the purity of the state of the remaining 
qubits in any given run (trial) of the algorithm. 

The quantum data-base search algorithm initially operates in an N = 2™-dimensional 
Hilbert space, where the basis vectors are represented as n-bit binary strings, \x n ) = 
\b x ,i, b x ,2 ■ ■ ■ b Xi7l ), where every b X)m takes two values and 1. The aim of the Grover search 
algorithm is to map an equal superposition state of all N basis vectors into a definite target 
object \x n ) = \h,b 2 - ■ -b n ). 

The search process on the n-qubit system is carried out by repeated action of a unitary 
operator G n which performs a change of sign of the target state amplitude and an inversion of 
all basis state amplitudes about their mean value [T|. The symmetry of the Grover operator 
G implies that the dynamics is restricted to a 2D Hilbert space spanned by the target state 
\x n ) and the fully symmetric state ^== J2 Xn \ x n)- As the target state is one of the terms in 
the sum over all basis states, it is convenient to span the 2D space by the orthogonal basis 
formed by the target state \x n ) and \s n ) = Ed„^i„ \ x n)- In this basis, G n amounts to 
a rotation by an angle 9 n = 2 sin~ 1 (l/\/ r N), 

G n = cos 9 n I — i sin 9 n a y , (1) 

where a y and I are the Pauli-?/ matrix and the identity operator respectively. If the ini- 
tial state p n (0) of the n-qubit register can be described within the 2D subspace, then k 
applications of G n yields the state 

p n (k)=G k n p n (0)(G^. (2) 

A complete transition from the symmetric state (which is readily prepared as a product 
state with every register qubit in the even superpositions state (|0) + |l))/\/2) to the target 
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state is obtained by |v N consecutive rotations, i.e., repeated operations of G n . This leads 
to the square root speed-up compared to a classical search through all N elements. 

A. Losses 

We now consider the situation where a number of qubits are lost after k applications of the 
Grover operator. If m qubits remain after the losses, the subsequent evolution takes place 
within a Hilbert space of dimension M = 2 m . Using the fact that both the target state and 
the fully symmetric state factorize, \x n ) = \x m ) <g> \x„- m ), ^= J2 X „ \ x ™) = J2 Xm \ x m) ® 
-^j^ 52x n - m \ x n-m), we can rewrite the basis vector \s n ) = [J2 Xn \x n ) - \x n )] in terms 

of the corresponding basis vectors of the subspaces containing m and n — m qubits, 

\Sn) fn,m\Sm) ® \Sn— m) 9n,m\%m) ® \Sn— m) ^n,m|^m) ® |£n— m) i (3) 

where / n , m = \j {N ~^ M ^ ~, 9n,m = sj m(n-i) and hn > m = \j ^=T- lt thus follows > that after 
discarding n — m qubits, the reduced density matrix of the remaining m qubits is given by 

Pm(k) = Tr n ^ m [p n {k)\ = (x n -rn\Pn(k)\x n -m) + ( s n-m \pn (k) |<S n _ m ) . (4) 

The key observation is that even though the system loses qubits, its effective dynamics can 
still be described by a density matrix evolution in the two-dimensional space spanned by 
the basis vectors \x m ) and \s m ). 

After performing the partial trace one finds that the elements p^{k) = (x m \p m (k)\x m ) , 
p£(k) = (s m \Pm(k)\x m ) and p%(k) = (s m \p m (k)\s m ) of the reduced density matrix p m (k) 
can be expressed in terms of the elements of the density matrix p n (k) as follows: 

= Pn(k) + gl mP f(k), f%(k) = (1 - gljpfik), (5) 

Pm( k ) = h n , m pn{k) + fn,m9n,mP^ (k). 

Since the lost qubits are entangled with the register one always finds, Tr[p^(/c)] < Tr[p^(k)], 
i.e., the reduced density matrix of the surviving qubits is always mixed. 

III. MASTER EQUATION 



As the system loses the qubits one by one in a random manner, the actual number of 
remaining qubits at any given time becomes a random variable, and our knowledge about the 
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system is appropriately described by a weighted sum over corresponding two-dimensional 
density matrices. We take the duration of a full Grover step as unit of time, and since each 
Grover step implements a small angular rotation (Eq. (pj)), we treat this as an effective 
continuous rotation in each m-qubit subspace, driven by the Hamiltonian H m = u m <j y , 
where u m = 2/a/M = 2~( m ~ 2 )/ 2 accounts for the rotation angle per time step in the space 
of m-qubits, and we set h = 1. The losses constitute a continuous random process with an 
average feeding of probability towards states with smaller qubit numbers. 

Let 7 denote the (small) probability for any qubit to be lost during the duration of a single 
Grover step. Due to the small values of the Grover rotation angle and the loss probability, 
we will represent both processes by a continuous time master equation 
dp (ij) 

— 7^ — = -i[H m , p m (t)\ - m>yp m (t) + (m + l)~/£(p m+1 (t)) (6) 

where p m {t) represents the projection of the density matrix on the subspace with m surviving 
qubits, and (m + 1)7 £(p m+1 (t)) represents the map in Eq. ^ for the loss of a single qubit, 
n — m + 1. The coupling due to losses of different projections of the density matrix is 
illustrated in Fig lb. 

The effect of the loss operator coupling the m th and the m + l st subspace density matrices 
can be found explicitly using the above equations, 
dp m (t) 



dt 
dt 



-m>yp m (t) + (m + l)>yp m+l (t) (7) 

2u m p u m - m 1P Z(t) + (m + 1) 7 ((1 - A m )pZ +l {t)) + A mPm+1 {t)) 

= -a^£ - m7p« $ + (ro + 

where, p m (t) = Tr[p m (t)], p^{t) = Pm(t) — p^(t) gives the population difference between 
the target and the symmetric states \x m ), \s m ), and p^if) = (p^(t) + p^(t))/2, gives the 
real part of their coherence in each subspace. The arguments leading to Eq. (|5) yield, for 
the loss of a single qubit, the coefficients A m = 2j J +1 and B m = ^/ 2 m"-i m 

Note that we used the same symbol for the density matrix in Eqs. (|4j [5]) which was 
derived for a definite number of qubits and therefore has unit trace, while the first equation 
111 Eq. for the weight distribution p m (t) leads, irrespective of the Grover dynamics, to 
the binomial distribution for the surviving number of qubits, 

p m (t)=( n )e- m *(l-e-*)™ (8) 



We now address the question of how much population has been transferred to the target 
state within the different subspaces, since this will tell us with what probability one obtains 
correct outcomes upon readout of the qubits remaining at the end of the calculation. 

Since there are only 2n coupled equations, one can solve the dynamics exactly by numer- 
ical integration for even very large numbers of qubits. The parameter regime of interest is 
one, where we may loose a significant fraction of the qubits, but of course not all of them, 
during execution of the Grover algorithm. Such a case is illustrated in Fig. 2, where we start 
with 24 qubits, and a loss rate leaving only 7-8 qubits at the anticipated end of the Grover 
search after ~ 3217 Grover operations. The figure shows that after about ~ 1300 Grover 
steps, half of the qubits are still present, and the different two-dimensional subspaces are 
approximately populating the fully mixed states, with equal weight on the target state and 
on the uniform state with the corresponding number of qubits, 

n 1 

P{t) ~ S ^2Prn{t)- [\x m )(x m \ + |s TO )(s m |] . (9) 
m=l 

Reading out the register at this time, thus provides on average ~ n/2 bit values, which are 
all correct with a nearly 50% probability (and which attain random and uncorrelated values 
with the remaining probability). 

Note that the number of surviving qubits is governed by the decay rate 7, while the 
approach towards the even mixture of the target and symmetric states of the surviving 
qubits is governed, in parts, by the m-dependent co m and the relaxation terms in Eq. <Q 
due to the trace over the lost qubits. 

To estimate this latter dynamics, we note that the binomial distribution is narrow with 
respect to the m-dependence of the coefficients in Eq. ([7]), and hence, we may evaluate the 
average dynamics of the m-qubit states by considering 

P"(*)=X>m(*)' P"(t)=XX(f), (10) 
m m 

and evaluating the different m-dependent coefficients in Eq. ^ at the time dependent 
average value 

m{t) = mp m (t) = ne~" yt . (11) 

m 

Summing over all m-subspaces and assuming M 3> 1, such that one can approximate the 
coefficients, A m w and B m « l/y/2 in Eq. one obtains the average dynamics, governed 
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FIG. 2: The relative probability T m (t) = {x m \p m \x m ) /p m (t) of the target state component in 
subspaces with 12 (red, steepest rise), 14 (green, second steepest rise) and 19 (blue, third steepest 
rise at small Nq) remaining qubits (out of n = 24), plotted as functions of the number of Grover 
steps Nq. The brown, dashed line shows the weighted probability F{t) = ^2 m PmJ^mit) that all 
bits match the target bit values. The results are shown for a loss rate of 7 ~ 4 x 10 -4 (per Grover 
step) leading to the average fraction of surviving qubits shown with the black line, while the target 
state probability in the ideal lossless case is shown with the gray line, reaching unity for Nq ~ 3217. 



by the two coupled equations 

dp w (t) 

dt 
dp u {t) 

dt 



2u(t)p u 

-2u(t)p w -T(t)p u . (12) 



These are the coupled equations for a two-level system with a time dependent Rabi-flopping 
rate, u(t) = 2-(™W _2 >/ 2 and dephasing rate T(t) = [(m(t) + l)B m{t) - fh(t)]j ~ ^ + 
(l - 3.) n^. 

Clearly if the Rabi-flopping is allowed to act long enough to distribute the population 
evenly among the states, before all the qubits have been lost, the state of the remaining 
qubits is the fully mixed state, with a 50% probability that all the bit values will be read at 
the correct target values. 
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IV. RECONSTRUCTION OF THE TARGET STATE 



From the above analysis it can be seen that we may loose a significant fraction of the 
qubits during the Grover search, and still, we may pick a moment in time, where, e.g., half 
of the qubits are left, and where, there is a 50% chance that their states all match the 
corresponding bit values in the correct target state. The probability that any single qubit 
is measured in its correct state is therefore 75%, and by repeating the protocol a number of 
times, so that we accumulate multiple read out values for every bit, a simple majority vote 
will yield the target bit values. 

We can, however, improve the statistical certainty of the result by exploiting the fact 
that with 50% probability, all the measured bit values are correct, while with the remaining 
50% probability they are all random. This implies that an m bit string yields m correct bits 
with a probability of 0.5 + 0.5 • ^ — > 0.5 for large m. Rather than estimating the correct bit 
values independently by majority votes, we suggest to use the correlations in the data and 
assign more weight to the data strings which agree with the majority votes in the largest 
number of bit positions and which are hence most likely to represent correct rather than 
random bits. 

We assume that we have obtained an experiment consisting of K different trials, where 
rrik qubits have survived in each trial, k — 1 • • • K. We represent the measurement outcomes 
of each trial as a string with n entries (b\ , • • • ,b^) attaining the values and 1 for the 
readout of a qubit in state |0) and |1), and for convenience the detection of a missing qubit 
is represented by the value 0.5 (see the inset in Fig. 3). 

For two output strings and (b\, ■ ■ ■ ,b*) and (b{, ■ ■ ■ ,b J n ), we perform a bit-wise compari- 
son, and introduce the quantifier 



otherwise 

v 

For each trial, we then define a measure for how well all its bit values agree with the 
values obtained in the all other (j ^ k) trials 



Xkj(i) 



S -1 



1 



(6*, 6*') = (0,0) or (1,1) 
(6^) = (1,0) or (0,1) 



n 




(13) 



j^k i=i 
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The values of Ck fall in two groups: random bits lead to values of Ck fluctuating around 0, 
while the positive contributions to Ck for correct bit strings from other correct strings cause 
a bias of their distribution towards positive values. Using this, we define a new, weighted 
estimate for the individual bit-values, 



(14) 

k ' ' 

The procedure does not suggest a preferred value for bit values when the argument of the 
sgn function vanishes, but this will occur with smaller and smaller probability the larger 
the number of trials. 

We have simulated trial records with half the qubits remaining on average, and after 
adequate time to render the density matrix close to the fully mixed state in each relevant 
two-dimensional m component, such that each string is with almost equal probability a 
random string and a correct string (in the places of remaining qubits). We have then 
analyzed numerically the outcome of the filtering algorithm presented, and we show results 
in Fig. 3. 

The calculations are done for a system with n = 24 qubits, and the blue curve (circles) 
shows how many bit values are correctly estimated by a simple majority vote after K = 10 
trials. This value fluctuates from experiment to experiment, and we show its probability 
distribution, according to which on average about 83% of the bits are correct and with 
about 10% probability all bit values are correct. The red curve (squares) shows the more 



advanced filtering, where the trial outcomes are assigned weights as described in Eq. (14). 
We observe a significant improvement in the average number of correct bits (95%) and with 
a 50% probability all the bit values are now determined correctly. Since the tasks solved by 
the Grover algorithm often entail a simple method to verify the outcome, the latter scheme 
presents a very efficient way to obtain the solution of the search problem in the presence of 
losses. 

In the numerical example we lost half of the qubits on average and we included K — 10 
trials in our analysis. With more trials we further reduce the probability that bit values 
remain unsettled but as indeed, a fraction of the trials predict correct bit values, only a 
finite number of trials is needed, and the overhead in performing the Grover algorithm that 
number of times to reconstruct the missing bits, does not change the y/N speed-up of the 
method. 
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FIG. 3: The probability distribution of the number n c of correctly estimated qubits for an initial 
register size of n = 24 qubits. In a given experiment we make 10 trials, to obtain the final bit 
string, and we compare the cases where the bit values are obtained (i) by a simple majority vote 



on the raw data (blue circles) (ii) by the correlated data approach (red squares) using Eq. (14). To 
obtain the probability distributions shown, we have averaged the outcome of 10 4 such experiments. 
Inset: The upper left corner of the measurement table displaying the results for each trial. The 
dark, green boxes indicate the lost qubits, represented by the numerical values b\ = 0.5 in Eqs. 



(13, 14). 



V. PHYSICAL IMPLEMENTATION WITH RYDBERG BLOCKADE 



It is a key assumptions in the above analysis, that the conditional phase shift operations 
that make the operator G, Eq. Q, work correctly within the subspace of the qubits actually 
present. This may not be true for a general implementation of quantum computing, where 
unitary operations on a register are accomplished through sequences of one- and two-bit 
operations, and where missing physical particles introduce errors that cascade into the time 
evolution of the ones being still present. 

In a recent quantum computing proposal [7J, it was shown that the Rydberg blockade 
interaction [131 [HI [T9142T] if it extends to any pair of atoms trapped in separate microtraps 
within a small volume in space, can be used to efficiently implement the Grover algorithm. 
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The basic ideas of the proposal is to use 7r-pulses to sequentially excite one atom after the 
other from a definite ground state or ground state superposition to the Rydberg state. As 
long as no atom occupies the ground state addressed by the laser pulse no excitation is 
formed, while the first atom being excited by a pulse will block any further excitation in 
the sample. Applying sequential excitation 7r-pulses followed by sequential de-excitation 
7r-pulses in the reverse order to the atoms, will return all population to the initial ground 
state, but a change of sign applies if a single atom was excited. This occurs if any non-zero 
number of atoms populates a state coupled by the laser field. If all excitation pulses are 
implemented from the qubit state different from the target bit value, this change of sign is 
readily encoding the Grover conditional phase. I.e., the target state component experiences 
a change of sign with respect to all other components having one or more erroneous qubits. 
The sign of the fully symmetric state with respect to all other states is equivalently obtained 
by 7r-pulse excitation and de-excitation sequences from the "erroneous" state (|0) — |l))/\/2~, 
leaving only the phase of the symmetric state unchanged. For details, see [7]. 

Missing atoms cannot be excited, and since the above operations work by putting a 
change of sign on any state, where one or more atoms can be excited from a wrong ground 
state or ground state superposition, the occasional missing of an atom does not disturb the 
execution of the algorithm. One may apply the laser pulse sequence as if all atoms are 
present, and the pulses illuminating empty traps have no effect on the implementation of 
the desired operations on the remaining atoms. 

VI. CONCLUSION 

In conclusion we have presented here an analysis of Grover search algorithm in the pres- 
ence of qubit losses. The effective dynamics in the presence of losses takes place in n 
two-dimensional subspaces, which allows a simple analytical and numerical investigation of 
its performance in the presence of losses. We showed that if up to half of the qubits are lost 
on the time scale of execution of the ideal Grover search, the resulting state of the remaining 
qubits has a large overlap with the desired target state restricted to the same qubits. We 
established a correlated data analysis method that performs better than a simple majority 
vote over repeated experiment, and we demonstrated that the full target state can be recon- 
structed by repeating the experiment a number of times that scales only linearly with the 
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number of bits in the register. I.e., the y/N quantum speed-up is not deteriorated by the 
losses. 

In schemes where it is possible to register when particles are lost, one may adapt any 
implementation of the Grover algorithm, as long as one adjusts the operations on the system 
such that the algorithmic steps are always applied correctly to the register composed of the 
qubits still present. 

We presented an experimental implementation scheme with neutral atoms and Rydberg 
blockade interactions, which naturally lends itself to demonstration of our results, because 
it can be applied on a lossy physical system without any attention to the losses during the 
operation of the algorithm. Current experiments with trapped arrays of atoms subject to 
Rydberg excitation J2DH22], may thus implement the protocols discussed in this paper and 
carry out basic tests of the Grover algorithm on multiple qubits, even beyond the time scale 
of individual particle losses. 
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